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A discrete-event approach, which has already been shown to give a cause-and-effect explanation of many 
quantum optics experiments, is applied to single-neutron interferometry experiments. The simulation algorithm 
yields a logically consistent description in terms of individual neutrons and does not require the knowledge of 
the solution of a wave equation. It is shown that the simulation method reproduces the results of several single- 
neutron interferometry experiments, including experiments which, in quantum theoretical language, involve 
entanglement. Our results demonstrate that classical (non-Hamiltonian) systems can exhibit correlations which 
in quantum theory are associated with interference and entanglement, also when all particles emitted by the 
source are accounted for. 
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I. INTRODUCTION 



Quantum theory has proven extraordinarily powerful for describing the statistical properties of a vast number of laboratory 
experiments. Conceptually, it is straightforward to use the quantum theoretical formalism to calculate numbers that can be 
compared with experimental data, at least if these numbers refer to statistical averages. However, a fundamental problem 
appears if an experiment provides access to the individual events that contribute to the statistical average. Prime examples are 
the single-electron two-slit experiment JJl], neutron interferometry experiments \%\ and similar experiments in optics where the 
click of the detector is identified with the arrival of a single photon yl. 

Although quantum theory provides a recipe to compute the frequencies for observing events it does not account for the 
observation of the individual detection events themselves flU)]. From the viewpoint of quantum theory, the central issue is how 
it can be that experiments yield definite answers. As stated by Leggett |@]: "In the final analysis, physics cannot forever refuse 
to give an account of how it is that we obtain definite results whenever we do a particular measurement". For a recent review 
of various approaches to the quantum measurement problem and an explanation of it within the statistical interpretation, see 
Ref.01. 

Perhaps the most simple and clear demonstration of the fundamental nature of this problem is provided by two-path interfer- 
ence experiments with electrons, photons, or neutrons. According to Feynman, the observation that the interference patterns are 
built up event-by-event is "impossible, absolutely impossible to explain in any classical way and has in it the heart of quantum 
mechanics. In reality it is the only mystery." ]8[]. 

Reading "any classical way" as "any classical Hamiltonian mechanics way", Feynman's statement maybe difficult to dispute. 
However, taking a broader view by allowing for dynamical systems that are outside the realm of classical Hamiltonian dynamics, 
it may be possible to model the gradual appearance of interference patterns through a discrete-event simulation that does not 
make reference to wave theory. This is precisely the approach taken in the present paper which is not about interpretations 
or extensions of quantum theory (see Ref. [2] for an overview) but adopts a new paradigm l9T Jl"lll to deal with the fact that 
experiments yield definite results. 

Feynman's statement that the event-by-event realization of an interference pattern is the only mystery suggests that creating 
interference patterns by a (local and causal) discrete-event process may be an important step in demystifying this aspect of 
quantum phenomena. Neutron interferometry is a close-to-ideal experimental technique to address this issue |2|]. The basic 
device used in the neutron interferometry experiments which are covered in this paper is a Laue-type interferometer B [l2l [l3ll . 
A large, perfect crystal of silicon is cut as shown in Fig. Q] The crystal plate BSO acts as a beam splitter: neutrons incident 
from the left are transmitted with or without being refracted by this plate. Neutrons refracted by beam splitters BS 1 and BS2 are 
directed to the third plate (BS3) which also acts as a beam splitter. Neutrons which are not refracted by beam splitters BS1 and 
BS2 leave the interferometer. To observe interference, the crystal planes of the different components have to be parallel to high 
accuracy lfl2ll and the whole device needs to be protected from vibrations ITill . All beam splitters are assumed to have the same 
reflection and transmission coefficients |2[]. Neutron detectors can have a very high, almost 100%, efficiency yfl. 

Many neutron interferometry experiments show that the intensity in the O- and H-beam, obtained by counting individual 
neutrons for a certain amount of time, exhibit sinusoidal variations as a function of the phase shift %, a prime characteristic 
of interference 0], Feynman's "mystery" pops up immediately if one wants to entertain the idea that only waves can produce 
interference. 

Adopting a wave-packet picture for an individual neutron, the wave packet first splits in two parts at BSO, then each part 
splits in two at BS1 and BS2. Two of the four parts go off to infinity, the other two parts "reunite" at BS3. At BS3 the 
merged wave packet splits again in two parts. Only one of these parts triggers a detector. It is indeed a mystery how four wave 
packets can conspire to do such things. Assuming that only a neutron, not merely a part of it can trigger the nuclear reaction 
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FIG. 1: Picture of the perfect crystal neutron interferometer 1 12]. BS0,...,BS3: beam splitters; phase shifter: aluminum foil; neutrons that are 
transmitted by BS1 or BS2 leave the interferometer and do not contribute to the interference signal. Detectors count the number of neutrons in 
the O- and H-beam. 



that causes the detector to "click", on elementary logical grounds, the argument that was just given rules out a wave-packet 
picture for the individual neutron (invoking the wave function collapse only adds to the mystery) but there is no conflict with the 
statistical interpretation of quantum mechanics JH0]. As long as we consider descriptions of the statistics of the experiment with 
many neutrons, we may still think of one single "probability" wave propagating through the interferometer and as the statistical 
interpretation of quantum theory is silent about single events, there is no conflict with logic either ||2[]. In this paper, we do not 
solve the aforementioned mystery but give an affirmative answer to the question whether it is possible to construct a logically 
consistent, cause-and-effect description in terms of discrete-event, particle-like processes which produce results that agree with 
those of neutron interferometry experiments and the quantum theory thereof. 

In previous work J^, [Til, fo - l2lll we have demonstrated, using an event-based corpuscular model, that interference is not 
necessarily a signature of the presence of waves of some kind but can also appear as the collective result of particles which at 
any time do not directly interact with each other. In general, the event-based approach takes as a starting point the observation 
that experiments yield definite results, such as for example the individual detector clicks that build up an interference pattern. 
We call these definite results "events". Instead of trying to fit the existence of these events in some formal, mathematical theory, 
in the event-based approach the paradigm is changed by directly searching for the rules that transform events into other events 
and, which by repeated application, yield frequency distributions of events that agree with those predicted by classical wave 
or quantum theory. Obviously, such rules cannot be derived from quantum theory or, as a matter of fact, of any theory that is 
probabilistic in nature simply because these theories do not entail a procedure (= algorithm) to produce events themselves. 

The paper is structured as follows. In Section [HI we specify the event-based model in detail. Sections IHI14VHI present 
our results for the basic neutron interferometer (see Fig. [TJ, experiments with stochastic and deterministic absorption, a Bell 
inequality test, an experiment that creates entanglement between the neutron path, spin and energy, and experiments that are 
performed in a non-stationary regime. For reference and to facilitate comparison, for each of the experiments that we discuss 
in this paper, we give the results of the quantum theoretical description of these experiments, adopting the terminology that is 
commonly used in quantum theory. In contrast, when we discuss the event-by-event, particle-like models of these experiments, 
there is no need to invoke concepts such as probability amplitudes, particle-wave duality etc. Our conclusions and outlook are 
given in Section rVIIII 

II. EVENT-BASED MODEL 

The event-based approach has successfully been used to perform discrete-event simulations of the single beam splitter and 
Mach-Zehnder interferometer ex peri ment of Grangier et al. ll22ll (see Refs. J^dUdll]), Wheeler's delayed choice experiment of 
Jacques et al. If23tl (see Refs. Il llll6l[l7ll ). the quantum eraser experiment of Schwindt et al. |24ll (see Ref. UHO), double-slit 
and two-beam single-photon interference experiments and the single-photon interference experiment with a Fresnel biprism of 
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Jacques etaJ. ill (see Ref. lfTT]JT 



Agafonov et al. B27I1 (see Ref. 
experiments of Aspect et 



LefTTJ 

aJ. E, 



]), quantum cryptography protocols (see Ref. 1261]). t he Hanbury Brown-Twiss experiment of 



, 20]), universal quantum computation (see Ref. [ 10, 28]), Einstein-Podolsky-Rosen-Bohm-type 
3C ] and of Weihs et al. i3~ltl (see Refs. ITTl |32 l- l37l p . and the propagation of electromagnetic 
plane waves through homogeneous thin films and stratified media (see Ref. fill 13811 ). An extensive review of the simulation 
method and its applications is given in Ref. Ill 111. Proposals for single-particle experiments to test specific aspects of the event- 
based approach are discussed in Refs. lfl9l l2lll . For many different optics experiments, the event-based corpuscular model 
reproduces the probability distributions of quantum theory or results of Maxwell's wave theory by assuming that photons have a 
particle character only. 

The event-based corpuscular model is free of paradoxes that result from the assumption that photons exhibit a dual, wave- 
particle behavior, and as we demonstrate in this paper, the same holds for neutrons as well. A crucial property of the event-based 
corpuscular models is that they reproduce various "wave results" observed in different experiments without any change to 
algorithms modeling the particles and components (e.g. beam splitters) lUlll . Although the event-based algorithms can be given 
an interpretation of a realistic cause-and-effect description that is free of logical difficulties, in the present stage of development 
it is difficult to decide whether or not such algorithms or modifications of them are realized by Nature. Only new, dedicated 
experiments may teach us more about this intriguing question. 



A. Definition of messenger and message 

A neutron is regarded as a messenger, carrying a message. As in our earlier event-based models for quantum optics experi- 
ments [ 1 1], we represent a message by the two-dimensional complex-valued unit vector 



e'> (I, cos(0/2) 
e'> (2) sin(0/2) 



(1) 



As is often the case, it is convenient though by no means essential to work with complex-valued vectors. The message Eq. ([]]) 
encodes the time of flight and the magnetic moment of the neutron. 

In a pictorial manner, the neutron carries with it a clock, the hand of which rotates with angular frequency v (to be discussed 
later). The clock may be used by event-based processors, mimicking the interaction of neutrons with materials, to determine 
the neutron's time of flight. Similarly, if we think of the neutron as a tiny classical magnet spinning around the direction 
m = (cos (j) sin 9 , sin (j) sin 9 , cos 9 ), relative to a fixed frame of reference defined by a magnetic field, then, the two angles (j) and 
9 suffice to specify the magnetic moment. 

According to Eq. ([]]), within the present model, it is postulated that the internal state of the neutron is fully determined by the 
three angles \jf^\ yf^ 2 \ and 9 and by rules, to be specified, by which these angles change as the neutron moves through space. In 
Eq. (Q~|), we have introduced three angles to characterize the message. The difference <j) = y/W — y/ 2 ' and 9 suffice to represent 
the magnetic moment and the third degree of freedom is used to account for the time of flight of the neutron. 

At this stage of the development, it is not clear whether the model of the messenger that we describe here is sufficient to explain 
all possible neutron interferometry experiments that might be carried out but to explain the neutron interferometry experiments 
which are covered in this paper, it cannot be simplified further. 

As the messenger moves for a time T, it is postulated that the message changes according to the rule 

y <- e ivT y, (2) 

where T is the time of flight, relative to the time of creation of the messenger, and v is an angular frequency. A monochro- 
matic beam of incident neutrons is assumed to consist of neutrons that all have the same value of v [2]. Put differently, of all 
the neutrons created in the fission process, the purpose of the monochromator is to select those neutrons that share the same 
characteristics: the velocity and direction in a classical mechanical picture [2] and the angular frequency v and direction in the 
event-based picture. The direction merely serves to send the selected neutrons to the interferometer. 

In the event-based picture, messengers can travel along a single path only. As they travel through the interferometer (one at a 
time) and are detected by one of the detectors, their times of flight may be different from messenger to messenger, depending on 
which path they followed and the delay they experienced in the material that acts as a phase shifter. Still within the event-based 
picture, the experimental fact that the measured intensity depends on the position of the phase shifter is a direct proof that the 
messenger conveys its time of flight to the processors. Hence it must have some kind of internal clock. 

A very plausible choice would be to relate v to the energy E of the neutron, that is we could make the hypothesis that 
v oc E/h where Planck's constant appears as a scale factor to render vT dimensionless. However, in this paper, the emphasis 
is on demonstrating that a particle-only model can reproduce the interference phenomena observed in neutron interferometry 
experiments and to simplify matters, we only consider idealized experiments with monochromatic beams of neutrons. In this 
case, the actual value of v does not affect the detector counts. Event-based simulations of experiments in which the actual 
value(s) of v are important, e.g. experiments which involve gravitation I2I I39I - I41I1 . are left for future research. 
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DLM 




FIG. 2: Diagram of a DLM-based processing unit that performs an event-based simulation of the beam splitters in the neutron interferometer 
(see Fig. [TJ. The processing unit consists of three stages: an input stage (DLM), a transformation stage and an output stage. The solid lines 
represent the input and output ports of the device. The presence of a message y (incident neutron) is indicated by an arrow on the corresponding 
port line. The DLM has storage for two real numbers (x) and two complex vectors Yo and Yj that are updated according to the rules Eqs. (|5j 
and I0, respectively. This data is combined to yield a 4-dimensional complex-valued vector which, after transformation by a matrix T, is fed 
into the output stage which decides through which port the (modified) message z leaves the device. The dashed lines indicate the data flow 
within the unit. 



In the presence of a magnetic field, a magnetic moment rotates about the direction of the magnetic field according to the 
standard, classical equation of motion. In terms of the message, this corresponds to a rotation of y about the same direction. 
As Eq. ([TJ suggests, the magnetic moment is represented through the well-known Bloch-sphere representation of a spin- 1/2 
particle ]5[]. Exploiting the relation between rotations in three-dimensional space and rotations in spin- 1/2 Hilbert space, in the 
presence of a magnetic field, the message changes according to the rule 

y <- i Xo*B*a>B,+o'Bj yi (3) 

where a x , O , and a z are the Pauli spin-matrices and B = (B x ,B y ,B z ) denotes the magnetic field vector. Although Eq. (01 is 
reminiscent of the rotation operator of a spin- 1/2 quantum object, in the present context Eq. Q is just a convenient construct to 
implement rotations in three-dimensional space. 



B. Particle source 



The source creates messengers and initializes the message. In order to demonstrate that the class of models which we consider 
can produce interference without solving wave equations, we explicitly exclude the possibility that at any time there is more than 
one messenger passing through the interferometer, an assumption which is often made in the discussion of neutron interferometry 
experiments IH2I1 . In the simulation, it is trivial to realize this condition: except for the first particle, the source creates a new 
particle only after the previous particle has been detected. It is also straightforward to let the source produce particles with 
specific properties. For instance, a fully coherent spin-polarized beam is simulated by generating messengers with the message 
given by Eq. ([TJ where l/A 1 ), y^ 2 ', and 9 are the same for all messages. Throughout this paper, the total number of particles 
generated by the source is denoted by N. 



C. Beam splitter 

In Fig. [2] we show the diagram of the event-based processor that simulates the operation of a beam splitter. This processor has 
three stages. The input stage consists of a so-called deterministic learning machine (DLM) Jgl [ill EH] • This machine is capable 
of learning, on the basis of the individual events, about the relative frequencies of messengers arriving on ports and 1. In 
neutron interferometry experiments, it is assumed that at any time, at most one neutron passes through the interferometer fl fl^l . 
In the event-based approach, this assumption implies that the DLM receives a message on either input port or 1, never on both 
ports simultaneously. 
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FIG. 3: Simulation data of the internal variable xo as a function of the number of received input events n, generated by Eq. (0. Initially xq = 0. 
For n = 0, ...,999, the input is either v = (1,0) (v = (0, 1)) with probability 0.8 (0.2). For n = 1000, ... ,2000, the input is either v = (1,0) 
(v = (0, 1)) with probability 0.2 (0.8). The horizontal lines represent the probabilities 0.8 and 0.2 The data of xq is shown as markers connected 
by thin lines. The running average of xy is shown as a thick line. At n = 1001, the value of the running average is set equal to the current value 
of xq. Solid squares and solid line: y = 0.99; Open circles and dashed line: y = 0.5. The data is plotted for every 20 input events. 



The arrival of a messenger at port or 1 is represented by the vectors v = (1,0) or v = (0, 1), respectively. A DLM that 
is capable of performing the desired task has an internal vector x = (xo,Xi), where xq +x\ < 1 and Xk > for all k = 0, 1. In 
addition to the internal vector x, the DLM should have two sets of two registers Y& = (J^, 1,1^2) to store the last message y that 
arrived at port k. Thus, the DLM has storage for exactly 10 real numbers. 

Upon receiving a messenger at input port k, the DLM performs the following steps: it copies the elements of message y in its 
internal register Y* 

Y* <- y (4) 

while leaving Yi_£ unchanged, and updates its internal vector x according to 

x <- yx + (1 - y)v. (5) 

It is easy to see that xo +X\ < 1 at all times. Each time a messenger arrives at one of the input ports, the DLM updates the values 
of the internal vector x and overwrites the values in the registers Y^. Thus, the machine can only store data of two messengers, 
not of all of them. 

The parameter < y < 1 affects the number of events the machine needs to adapt to a new situation, that is when the ratio of 
particles on paths and 1 changes. By reducing y, the number of events needed to adapt decreases but the accuracy with which 
the machine reproduces the ratio also decreases. In the limit that y = 0, the machine learns nothing: it simply echoes the last 
message that it received iHdH]. If Y — > l - , the machine learns slowly and reproduces accurately the ratio of particles that enter 
via port and 1 . It is in this case that the machine can be used to reproduce, event-by-event, the interference patterns that are 
characteristic of quantum phenomena ll^ fTllfT5ll . 

For later applications, it may be useful to have some insight into the dynamics of this DLM. In Fig. [3] we show some 
representative results obtained by executing the rule Eq. (fSJ- From Fig.[3j it is clear that all the features that we have discussed 
are present in the data. As y = 0.99 is close to one, the processor learns slowly. It takes several hundreds of input events before xo 
fluctuates around the probability 0.8 for a (1,0) input event. If we change the latter from 0.8 to 0.2, the DLM reacts immediately 
but again it takes a few hundred steps to reach the stationary state that corresponds to the new input sequence. The relatively 
slow pace with which the DLM responds to a change of the input sequence has a significant impact on the running average, 
represented by the thick solid line. For the number of events shown, for y = 0.99, the running average does not come close to its 
asymptotic value (for 20000 input events it does, data not shown). For y = 0.5, we see that the DLM responds very fast but this 
at the cost of large fluctuations of xq. 
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For applications to the event-based simulation of quantum phenomena, the value of the running average is of no importance 
but the fluctuations of xq are. We will see later that large fluctuations reduce the visibility of the interference signal. The 
number of events, required to establish the stationary state, becomes an important issue for simulating experiments in which the 
conditions rapidly change with time (see Section lVHBI l. Otherwise it is not an issue. Summarizing: the parameter / determines 
the "quality" of the event-by-event model of the interferometer, the ideal interferometer corresponding to y— !> I - . 

Returning to the diagram of the processor, the second stage accepts a message from the input stage and transforms it into a 
new message. From the description of the DLM, it is clear that the internal registers Yo and Yi contain the last message that 
arrived on input port and 1 respectively. First, this data is combined with the data of the internal vector x, the components of 
which converge (after many events have been processed) to the relative frequencies with which the messengers arrive on port 
and 1, respectively. The output message generated by the transformation stage is 
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(6) 



where the reflection R and transmission T = 1 — R are real numbers that are considered to be parameters, to be determined 
from experiment. Note that in contrast to optics Hill where S- and P-polarized waves may behave differently upon reflec- 
tion/transmission [43], in the case of neutrons, the first matrix in Eq. (0 (reading from left to right) treats the first and second 
pair of the four-dimensional vector on equal footing, in concert with the quantum theoretical treatment in section El Al Further 
note that asxo+xi < 1 at all times and ||Y || = ||Yi|| = 1, we have |Z ,i | 2 + \Z . 2 \ 2 + |Z U | 2 + |Zi, 2 | 2 = 1. 

The output stage uses the data provided by the transformation stage to decide through which of the two ports a messenger 



(representing a neutron) will be sent. The rule is very simple. We compute z 
the rule 

k' = ®{z-M), 



*\2\ and select the output port k' by 



(7) 



where ©(.) is the unit step function and < ffl < 1 is a uniform pseudo-random number (which changes with each messenger 
processed). From a simulation point of view, there is nothing special about using pseudo-random numbers. On a digital computer, 
pseudo-random numbers are generated by deterministic processes and therefore the pseudo-random number generator may be 
replaced by any algorithm that selects the output port in a systematic, uniform manner fl fllll . as long as the zero's and one's 
occur with a ratio determined by z. In fact, we use pseudo-random numbers to mimic the apparent unpredictability of the 
experimental data only. 

The messenger leaves through either port k' — or port k' = 1 carrying the message 



^,1 



z kf,l 

Z k'.2 



(8) 



which, for internal consistency and modularity of the event-based approach, is also a unit vector. 



D. Detector 



In the simulation model, we simply count all neutrons that leave the apparatus through the O- and H-beam. In other words, 
we assume that the detectors have 100% detection efficiency. Note that real neutron detectors can have efficiencies of 99% and 
more Hll . 



E. What makes it work? 



Anticipating that the event-based processor described in this section will perform as expected, that is, produce the expected 
interference patterns, it may be useful to have a deeper understanding of how it can be that these patterns appear without solving 
a wave problem. 

Let us consider BS3 in Fig. [1] the beam splitter at which, in a wave picture, the two beams join to produce interference. 
The event-based processor simulating a beam splitter requires two pieces of information to send out particles such that their 
distribution matches the wave-mechanical description of the beam splitter. First, it needs an estimate of the ratio of particle 
currents in the O- and H-beam, respectively. Second, it needs to have information about the time of flight along the two different 
paths. 
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FIG. 4: Diagram of the interferometer shown in Fig.Q] BS0,...,BS3: beam splitters; <j>o and <j>i: phase shifters. Detectors count all particles 
that leave the interferometer via the O- and H-beam. In experiment and in the event-based simulation, neutrons enter the interferometer via 
the path labeled 1*0 only. The paths labeled *Pi, 1*2 > an d 1*3 are use d m the quantum theoretical treatment only (see text). Particles leaving the 
interferometer via the dotted lines are not counted. 



The first piece of information is provided for by the internal vector x. As explained above, through the update rule Eq. (0, 
for a stationary sequence of input events, x = (xq,x\) converges to the average of the number of events on input ports and 1, 
respectively. Thus, the intensities of the waves in the two input beams are encoded in the vector x. Note that this information is 
accurate only if the sequence of input events is stationary. 

After one neutron arrived at port and another one arrived at port 1, the second piece of information is always available in the 
registers Yo and Yi. This information plays the role of the phase of the waves in the two input beams. 

It is now easy to see that all the information (intensity and phase) is available to compute the probability for sending out 
particles according to the distribution that we know from wave mechanics. Indeed, in the stationary state, Eq. (0 is identical to 
the transformation of the wave amplitudes which we know from wave theory of a beam splitter 02l 14311 . 

The idea that the event-based model of a beam splitter has some memory and a learning capability may seem strange enough 
to reject the model at first sight. However, applying the same logic to for instance Maxwell's theory of electrodynamics, one 
should reject this model as well. Indeed, the interaction of the electromagnetic wave and a material invariably takes a form that 
involves memory. This can be seen as follows. In Maxwell's theory, for electromagnetic radiation with frequency CO, the (linear 
part of the) interaction of the electric field E(o) and a material takes the form P(co) = i~i(co)E(co) where P(co) and r\{co) are 
the polarization and dielectric susceptibility of the material, respectively B43I1 . Transforming this relation to the time domain and 
assuming that E(t = 0) = P(f = 0) = yields JH 

P(/) = / ri(t-u)E(u)du, (9) 
Jo 

where the memory kernel 77(f) is the Fourier transform of rj(co). Clearly, Eq. (|9]i shows that the response of the polarization 
vector to the electric field involves memory. 

It is instructive to make the analogy with the update rule Eq. (|5j more explicit. Assume that x^ and are the values of 
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time-dependent vectors x(f ) and v(f ) sampled at regular time intervals x. If x(f ) allows a Taylor series expansion, we may write 
Xk = x(xk), Xk-i — x(xk) — Tdx(t)/dt\ t=r i; + ff(x 2 ) such that the update rule Eq. (0 can be expressed as 

^ = -^- r x(0 + ^- 7 v(0. (10) 

dt xy xy 

In order that Eq. ( TTOb makes sense for x —> 0, we must have lim T ^o(l — 7)1^7 = r. This requirement is trivially satisfied by 
putting y = 1/(1 + xT). Then Eq. (TTOb takes the form of the first-order linear differential equation 

^ = -rx(0 + rv(r). (ii) 

Assuming x(0) = 0, the formal solution of Eq. (fTTT i reads 

x(f) = T f e-" r \(t-u)du 7 (12) 
Jo 

which has the same structure as Eq. (O. From the derivation of Eq. ( fTTI ). it follows that if we interpret x as the time interval 
between two successive messages and let x approach zero, then 7=1/(1+ xF) approaches one and the DLM defined by the 
update rule Eq. (0 "solves" the differential equation Eq. (fTTT i. Therefore, we may view Eq. (fTTT i as a course-grained, continuum 
approximation to the event-by-event process defined by Eq. (0. 

Summarizing, the general idea that objects retain some "memory" about their interaction with external agents (particles, 
fields,...) is not only common but even essential to some of the most successful theories of physical phenomena and can 
therefore not be used as an argument to dismiss a particular class of models. Furthermore, it is worth noting that Eq. (0 is not 
the only update rule which yields an event-based model that reproduces the averages predicted by quantum theory l9l. lllll . In 
other words, there is nothing "unique" to Eq. (0. Whether an event-based model accounts for what is actually happening on the 
level of single events can only be decided by experiments that address this specific question. 



III. NEUTRON INTERFEROMETER 



A detailed wave-mechanical description of the diffraction of neutrons by the perfect silicon plate and the complete neutron 
interferometer is given in Ref. J2|. In this paper, to simplify matters without giving in on the fundamental issues, we adopt an 
effective model for the scattering process of the neutron and the plate. We assume that the neutrons are monochromatic and 
satisfy the Bragg condition for scattering by the silicon plate [2]. This is not an essential simplification. In the theory of neutron 
interferometry, it is customary to compute the incoherent average over slight deviations from the exact Bragg condition and 
neutron energy [2] and the same can be done in the event-based approach as well (see Section lVHBI ). Thus, we will characterize 
the beam splitters BS0,...,BS3 by effective reflection and transmission coefficients r and f, respectively. 

Once it has been established that the event-based approach reproduces the results of wave theory, a ray-tracing scheme such 
as the one outlined in Ref. [2] can be combined with the event-based processors to yield a more complete description of how the 
individual neutrons propagate through the interferometer and produce interference. We leave this technically challenging topic 
for future research. 



A. Quantum theory 

A detailed quantum mechanical treatment of the interferometer depicted in Fig. Q] is given in Ref. M45I1 . Assuming that the 
incident wave satisfies the Bragg condition for scattering by the first crystal plate (BSO), the Laue-type interferometer acts as a 
two-path interferometer 14611 . The two-path interferometer may be represented by a more abstract, theoretical model, the diagram 
of which is shown in Fig. [4] This diagram is similar to the one of the Mach-Zehnder interferometer for light l43ll . except that the 
latter has mirrors instead of beam splitters BS1 and BS2. 

Quantum theory describes the statistics of the interferometry experiment depicted in Figs. [T]and|4]in terms of the state vector 

\V) = (^or^oi^n^n^^n^^if, (13) 

where the components of this vector represent the complex-valued amplitudes of the wave function. The first subscript labels 
the pathway and the second subscript denotes the direction of the magnetic moment relative to some B-field, the direction of 
which becomes relevant if the experimental outcome depends on the magnetic moment of the neutron (see later). This is not the 
case for the experiment shown in Fig. [1] hence there is no need to dwell on this aspect any further. As usual, the state vector 
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is assumed to be normalized, meaning that (^l 1 ?) = 1. Note that in the abstract representations of the experiments, such as in 
Fig. |4]for example, we use the notation *Pj — ) for j = 0, . . . , 3. 

As the state vector propagates through the interferometer, it changes according to 



r r \ / r r 
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m, (14) 



where t and r denote the transmission and reflection coefficients, respectively, and the subscripts i, j refer to the pair of elements 
of the eight-dimensional vector on which the matrix acts. Conservation of probability demands that \t\ 2 + \r\ 2 = 1. For a later 
application, we have included in Eq. (TPfl i. a path-dependent absorption parameterized by the coefficients a (0 < a < 1) and b 
(0<Z><1). 

In neutron interferometry experiments, particles enter the interferometer via the path corresponding to the amplitude on ly 
(see Fig.[U), meaning that \ x i l ) = (1,0,0,0,0,0,0,0). The probabilities to observe a particle leaving the interferometer in the H- 
and O-beam are then given by 

PH = \% t \ 2 + \% i \ 2 = R (aT 2 + bR 2 -2RTV^b cos x), (15) 

Po=|^ 3t | 2 + l^3|| 2 = R 2 T (a + b + 2V^bcosx) , (16) 

where % = <j>\ — <j>o is the relative phase shift, R — \r\ 2 and T = \t \ 2 = 1 — R. Note that pa and po do not depend on the imaginary 
part of t or r, leaving only one free model parameter (e.g. R). In the case of a 50-50 beam splitter (T = R = 1 /2) and zero 
absorption (a = b = 1), Eqs. (TT3T > and (fT6l i reduce to the familiar expressions pn = (l/2)sin 2 %j2 and po = (l/2)cos 2 j/2, 
respectively. The extra factor two is due to the fact that one half of all incoming neutrons, that is the neutrons that are transmitted 
by BS1 or BS2 (see Fig. [TJ, leave the interferometer without being counted. 

The expression Eq. (TTST i shows that the normalized (to the maximum value) O-beam intensity does not depend on the value of 
the reflection R. Furthermore, it follows from Eq. ( TToT l that the visibility of the O-beam is given by 



_ max xPo -mm xP o 2^/ab 

V{a,b) = ■ : = — — , (17) 

max^ po + mm x po a + b 



and that the modulation amplitude of the interference fringes is given by (a + b)V{a,b) /2 = vab 



B. Interferometer: Event-by-event simulation model 

Using the event-based processor described in Section|II] it is straightforward to construct a simulation model for the interfer- 
ometer shown in Fig. [4] Without any modification, we use the event-based model of a beam splitter to simulate the operation 
of BS0, BS1, BS2, and BS3. Neutrons that are not refracted by BS1 or BS2 leave the apparatus and do not contribute to the 
detection counts in the O- or H-beam. During their flight from BS1 or BS2 to BS3, the neutrons pass through a metal foil which 
changes their time of flight JUl. In the event-based model this effect of the metal foil is accounted for by the phase shifters 0o 
and 0i, see Fig. [4] Thereby it is assumed that the absorption of neutrons by the metal foil is negligible [2]. When the messenger 
passes through the phase shifter, its message changes according to 

y*-e»Jy, (18) 

where <j>j represents the change in the time of flight as the neutron passes through the metal foil on its way from BS1 to BS3 
( j = 0) or BS2 to BS3 ( j = 1 ). In neutron interferometry experiments, minute rotations of the foils about an axis perpendicular to 
the base plane of the interferometer induce large variations in <$>j 12j,|47[]. All the neutrons which emerge from the interferometer 
through the O- or H-beam contribute to the neutron count in these beams. 



C. Simulation results 



The simulation results presented in Fig.[5ja) demonstrate that the event-by-event simulation reproduces the results of quantum 
theory if 7 approaches one [Joj, [TTJ, [T5|] . Indeed, there is excellent agreement with quantum theory. In this example, the reflection 
coefficient of the beam splitters is taken to be R = 0.2. The parameter 7 which controls the learning pace of the DLM-based 
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FIG. 5: Event-by-event simulation results of the number of neutrons leaving the interferometer via the H-beam (red circles) and O-beam (blue 
squares) as a function of the phase difference % between the two paths inside the interferometer. For each value of x, the number of particles 
generated in the simulation is N = 100000. The lines are the predictions of quantum theory for a — b = 1. Solid line: pn, see Eq. d!5t : dotted 
line: po, see Eq. d!6t . (a) Model parameters: reflection R = 0.2, y = 0.99. (b) Same as (a) except that y = 0.5, reducing the accuracy and 
increasing the response time of the DLM. (c) Same as (a) except that to mimic the partial coherence of the incident neutron beam, the initial 
message carried by each particle has been modified by adding to l/A 1 ) and i|/ 2 ) a random angle drawn uniformly from the interval [—60°, 60°], 
reducing the amplitude of the interference, (d) Same as (c) except that y = 0.5. 

processor can be used to account for imperfections of the neutron interferometer. This is illustrated in Fig. |5jb) which shows 
simulation results for 7 = 0.5. 

The quantum theoretical treatment of Section ITlI Al assumes a fully coherent beam of neutrons. In the event-based approach, 
the case of a coherent beam may be simulated by assuming that the degree of freedom that accounts for the time of flight of 
the neutron takes the same initial value each time a message is created. In the event-based approach, we can mimic a partially 
coherent beam by simply adding some random noise to the message, that is when a message is created, a pseudo-random number 
in a specified range is added to y/ 1 ' and V/' 2 '. In Fig. EJc), we present simulation results for the case that the random angle is 
drawn randomly and uniformly from the interval \—n/3,n/3], showing that reducing the coherence of the beam reduces the 
visibility, as expected on the basis of wave theory 04311 . Comparing Fig. [2b) and Fig.|5jc), we conclude that the same reduced 
visibility can be obtained by either reducing 7 or by adding noise to the messages. On the basis of this interferometry experiment 
alone, it is difficult to exclusively attribute the cause of a reduced visibility to one of these mechanisms. For completeness, 
Fig. |5jd) shows the combined effect of decreasing 7 and adding noise to the messages on the visibility of the interference 
fringes. 

Conclusive evidence that the event-based model reproduces the results of a single-neutron interferometry experiment comes 
from comparing simulation data with experimental data. In Fig. [6] we present such a comparison using experimental data 
extracted from Fig. 2 of Ref. [ 14]. It was not necessary to try to make the best fit: the parameters R and 7 and the offset in % 
were varied by hand. As shown in Fig. [6] the event-based simulation model reproduces, quantitatively, the experimental results 
reported in Fig. 2 of Ref. H- 
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FIG. 6: Comparison between the data of a neutron interferometry experiment [14] (open symbols) and the results of an event-by-event 
simulation (solid symbols). Open circles: counts per second and per square cm in the O-beam; open squares: counts per second and per square 
cm in the H-beam; solid circles: number of particles per sample leaving the interferometer via path 0; solid squares: number of particles per 
sample leaving the interferometer via path 1. The experimental data has been extracted from Fig. 2 of Ref. 14. The simulation parameters 
R — 0.22 and y = 0.5 have been adjusted by hand to obtain a good fit and the number of incident particles in the simulation is N = 22727. 
Lines through the data points are guides to the eye. 



IV. STOCHASTIC AND DETERMINISTIC BEAM ATTENUATION 

The second series of experiments that we consider are neutron interferometry experiments in which the beam from BS0 to 
BS1 is attenuated either by a partial, stochastic absorber JHHgl], (see Fig.|7ta)), an absorbing lattice | | p,|49ll50ll or by a chopper, a 
rotating absorbing disc (see Fig.|7Jb)) that periodically blocks neutrons from traveling to BS1 ll^L 1481 - 45011 . We denote the average 
fraction of neutrons which pass the absorber/lattice/chopper by < a < 1 . 

Let us assume that the incident flux of neutrons is constant in time. Of all neutrons passing through the stochastic absorber, 
only the fractor a "survives" the interaction with the absorber material. Considering the case of the experiment with a chopper, as 
the rotation frequency of the chopper increases up to the point that each individual neutron passes the chopper with probability 
a, the difference between the deterministic and stochastic absorption is expected to disappear 0H45JJ]. Although in the 
experiment, the chopper rotates inside an aluminum chamber, the vibrations associated with the rotation have an adverse effect 
on the amplitude of the interference fringes and hence, the rotation frequency of the chopper was effectively limited to a few 
rotations per second [48]. This problem could partially be alleviated by using an absorbing lattice instead of a rotating disc [2, 
r49l l50ll . In this section, we assume that the rotation of the chopper does not affect the experimental outcomes other than by 
stopping particles from reaching BS 1 . 

A. Quantum theory 

The case of the stochastic absorber has already been treated in Section ITlIAl When the stochastic absorber is replaced by a 
chopper that either passes the neutrons or blocks them completely, quantum theory prescribes that the experiment is described 
in terms of a mixed state [5], that is the observed intensity is the weighted sum of the two different experimental situations. 
Denoting the average fraction of the neutrons which pass the chopper by a, the probability and visibility of the O-beam are given 
by [48]" 




p' = (1 — a)po(a = 0,b = 1) +apo(a = l,b = 1) = TR 2 (1 + a + 2acosx) 



(19) 
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FIG. 7: Diagram of the single-neutron interferometry experiment with a stochastic (left) and deterministic (right) absorber. 
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FIG. 8: Diagram of the single-neutron interferometry experiment with a stochastic or deterministic absorber 1 12]. BS0,...,BS3: beam splitters; 
00 and 0i : phase shifters. Neutrons pass through the stochastic or deterministic absorber with probability a. 



and 



V(a) 



max * Po- min xPo _ _ 
max x p' +min x p' 1 



2a 



(20) 



respectively. In contrast to the case of the stochastic absorber where the modulation amplitude of the interference fringes is given 
by i/a, when a chopper is used to block neutrons that travel towards BS1, the modulation amplitude of the interference fringes 
is given by (1 +a)V(a)/2 = a, that is it is linear in a 111 211 . 

In quantum theory, the difference between the stochastic and deterministic absorber enters through the choice of the state. In 
the former case, we use a pure state to describe the interference pattern whereas in the latter case, a mixed state which adds the 
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FIG. 9: Modulation amplitude of the intensity as a function of the probability a that the neutrons pass the absorber, obtained from event-by- 
event simulations of single-neutron interferometry experiments [48-50] with a stochastic (a) and a deterministic absorber (b,c,d) for different 
values of iV pc , the number of incident particles per cycle and per value of the phase shift %, and N c , the number of cycles, keeping the total 
number N pc N c fixed. Solid line: square root dependence predicted by quantum theory for the case of a pure state (see text); dashed line: 
linear dependence predicted by quantum theory for the case of a mixed state (see text); solid triangles: N pc = 1000, N c = 10; open squares: 
A'pc = 100, N c = 100; solid squares: N pc = 250, N c = 40; open circles: N pc = 500, N c = 20; solid circles: AT pc = 1000, N c = 10; open triangles: 
A'pc = 10000, N c = 1. Model parameters: reflection R = 0.2, y = 0.98 (a,b), y= 0.9 (c), y= 0.5 (d). Within the statistical fluctuations, 
the simulation results for the stochastic absorber do not depend on A'pc or N c . For all cases, increasing the number of cycles jV c reduces the 
statistical fluctuations only. The maximum count of neutrons in the O-beam is less than 2000. 



probabilities for the two different experimental configurations (chopper blocking or not). However, this stationary-state wave 
theory does not include as a parameter the rotation speed of the chopper. Only the fraction (a) of neutrons which pass with each 
cycle enters, not how many neutrons pass through one opening in the disc at a time. This problem calls for a solution of the 
time-dependent Schrodinger equation but, to the authors knowledge, there is no report of successful work in this direction. Put 
differently, in the quantum theoretical description of this experiment, the parameter (chopper rotation frequency) which connects 
Eq. (IT9b to Eq. ( fT6l > is lacking. 

Regarding the individual neutron as a particle, consider the case in which the chopper rotates very slowly relative to the pace 
with which the neutrons arrive. Then, when the position of the chopper allows the neutrons to pass, many of them pass before 
the chopper closes and during this period, we may expect to see the interference signal which is characteristic of the two-path 
interferometer. If the chopper blocks the beam, there is no interference. It is in this case that the mixed state describes the 
statistics of the experiment. 

But what if the chopper rotates very fast such that with each open/close change, on average one neutron impinges on the 
chopper? In this case, the detected signal is not simply the sum of two independent experiments (one with the path from BS0 
to BS 1 blocked and another one with no blocked paths) but should be the same as in the case of a stochastic absorber (with 
the corresponding value of a). The next subsection shows that the event-based model effortlessly reproduces this behavior and 
provides a unified, logically consistent description of these experiments. 
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B. Event-based model 

Both the case of a stochastic and time dependent absorber are easily incorporated in the event-based model of the interferome- 
ter. In the former case, particles leaving BSO towards BS 1 pass through the absorber if 2% < a where a is the fraction of particles 
that passes and M is a uniform pseudo-random number (which changes with each particle). In the latter case, the procedure is as 
follows. First, we define the (dimensionless) unit time interval by the time it takes for the chopper to open and close. The number 
of neutrons incident on the interferometer per unit time interval and the number of such intervals will be denoted by Ns and Nj, 
respectively. Thus, in one simulation run, the total number of neutrons created by the source is N = NsN]. Uniform pseudo- 
random numbers are used to generate N times in the interval [0,Afy]. These are the times, relative to the motion of the chopper, 
at which the neutrons will arrive at the chopper, if they followed the path from BSO to BS1. To a very good approximation, the 
interarrival times are distributed exponentially or, in other words, the events are created according to a Poisson distribution [52]. 
Neutrons are sent to the interferometer, one at a time, in chronological order. When a neutron arrives at the chopper, its arrival 
time is used to determine if the chopper is open or closed. If the chopper is open, the neutron continues its journey to BS1. 
Otherwise, it is removed from the system and does not contribute to the detection counts. All neutrons that appear in the O- or 
H-beam are recorded by the detectors. 

The data is collected as follows. First we choose < a < 1. Then, for each setting of the phase shift %, the source sends 
N = N VC N C neutrons to the interferometer. Here, Ap C is the number of incident particles per cycle and per value of the phase shift 
X and N c is the number of cycles. For a fair comparison between different cases, we keep N fixed. Increasing N c while keeping 
N pc constant reduces the statistical fluctuations only. As in the case of the interferometer without absorbers, the event-based 
processors in BSO, BS1, BS2 and BS3 perform their task, yielding counts in the O- and H-beam. By changing % an d repeating 
the simulation, we determine the visibility V(a). Plotting (1 + a)V(a)/2 as a function of a, we can directly compare to the 
experimental results [2, 48-50]. According to Eqs. ( TP71 ) and d20l . (1 +a)V(a)/2 is proportional to ^/a and a for a pure and 
mixed state, respectively. 

C. Simulation results 

In the event-based approach, a certain number of particles has to pass through the interferometer before an interference pattern 
appears (obviously, as in experiment, a single particle does not produce an interference pattern). In addition, with each change 
of the chopper position, the learning machine in BS3 has to adapt to the change in input as particles arrive on either port or 
1 when the path via BS 1 is not blocked and no particles arrive on port when the path via BS 1 is blocked. Thus, the learning 
machine should adapt quickly to a new situation, a requirement which is in conflict with the desire to reproduce the results of 
quantum theory, which demands y — > 1~. From these simple observations, we expect that the parameter yean be used, not only 
to control the visibility of the interference fringes (see Fig. |5ja,b)) but can also be used to produce features which in quantum 
theory, are characteristic of the mixed state. The results presented in Fig.|9]confirm these expectations. 

In Fig.|9ja), we present the results for the stochastic absorber. Disregarding the statistical fluctuations, the data nicely follow 
the y/a curve predicted by quantum theory and is in agreement with experiment 1I2I l48Tl50ll . Note that in this case, as long as 
A^ = Np C N c is fixed, the values of N vc and N c themselves should, and also do not matter because the process to block neutrons 
from reaching BS 1 is time independent. 

Simulation results for the case of the deterministic absorber are presented in Figs.|3b,c,d), corresponding to 7 = 0.98,0.9,0.5, 
respectively. For 7 = 0.98, see Fig. |9fb), the machines learn fairly slowly. Yet for N pc = 10000 (the value of N c merely affects 
the statistical fluctuations), the signature of the mixed state, the linear dependence of (1 +a)V(a)/2 on a, starts to appear 
(open triangles). For J — 0.9, see Fig. |9jc), the simulation produces the results of both the pure and the mixed state (solid and 
open triangles, respectively). Note that the event-based model also delivers intermediate results (solid circles), as observed in 
experiment HI |4§l For 7 = 0.5, see Fig. |9|d), the processors can quickly adapt to a new situation, yielding results that 
interpolate smoothly between the linear- and square-law dependence on a. 

Summarizing, the event-based model produces results that agree with the stochastic absorber or very fast chopper, a fast 
chopper and a slow chopper, but unlike in the quantum theoretical descriptions, without any modification to the simulation 
model but only by changing the number of neutrons per open/close cycle, that is the rotation frequency of the chopper relative 
to the number of incident neutrons. 



V. VIOLATION OF A BELL INEQUALITY 

The neutron interferometry experiment reported in Ref. 115 311 demonstrates that it is feasible to manipulate independently, the 
spatial and spin degree of freedom of massive particles. The experiment shows that it is possible to create correlations between 
these two degrees of freedom which, within quantum theory, cannot be described by a product state. The direct experimental 
evidence is that the data for this correlation violates a Bell-CHSH inequality [53]. 
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FIG. 10: Diagram of the single-neutron interferometry experiment to observe correlations that cannot be accounted for by a quantum system 
in the product state, see also Fig. 1 in Ref. 15311 . BS0,...,BS3: beam splitters. The combination of the mu-metal spin turner, phase shifter, and 
spin rotator allow the independent manipulation of the neutron's spatial and magnetic degrees of freedom. 



In this section, we show that the event-based model faithfully reproduces all the features of quantum theory for this exper- 
iment and, b y chang ing the model parameter y, can also reproduce the numerical values of the correlations, as measured in 
experiments 11531 [5411 . 

A diagram of the single-neutron interferometry experiment is shown in Fig. [10] see also Fig. 1 of Ref. lf53ll . Incident neutrons 
pass through a magnetic-prism polarizer (not shown) that produces two spatially separated beams of neutrons with their magnetic 
moments aligned parallel (spin up), respectively anti -parallel (spin down) with respect to the magnetic axis of the polarizer which 
is parallel to the guiding field B. The spin-up neutrons impinge on a silicon-perfect-crystal interferometer. On leaving beam 
splitter BSO, neutrons may or may not experience refraction. A mu-metal spin-turner changes the orientation of the magnetic 
moment from parallel to perpendicular to the guiding field B. In detail, the result of passing through the mu-metal spin-turner 
is that the magnetic moment of a neutron that travels towards BS1 (BS2) rotates by 7r/2 (-7C/2) about the y-axis. Before the 
different paths join at the entrance plane of beam splitter BS3, a difference between the times of flight (corresponding to a phase 
in the wave mechanical description) along the two paths can be manipulated by a phase shifter. The neutrons that experience 
two refraction events when passing through the interferometer form the O-beam and are analyzed by sending them through a 
spin rotator and a Heusler spin analyzer. If necessary, to induce an extra spin rotation of n, a spin flipper is placed between 
the interferometer and the spin rotator. The neutrons that are selected by the Heusler spin analyzer are counted with a neutron 
detector (not shown) that has a very high efficiency (« 99%) 15311 . 



A. Quantum theory 

An essential feature of the experiment is that as the neutron passes through the interferometer, its path and its magnetic moment 
become correlated ll53tl . The quantum theoretical description of the experiment reported in Ref. ll53tl requires a four-state system 
for the path and another two-state system to account for the spin- 1/2 degree-of-freedom. Thus, the statistics of the experimental 
data is described by the state vector Eq. (IT3T >. 

In the experiment [53], the neutrons that enter the interferometer have their spins up, relative to the direction of the guiding 
field B (see Fig. |T0]l. Thus, the state describing the incident neutrons is |\P) = (l,0,0,0,0,0,0,0) r , omitting irrelevant phase 
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FIG. 11: Quantum theoretical model of the single-neutron Bell-inequality test experiment 1 5 3] . Polarizer spin-up neutrons are injected in the 
interferometer. BS0,...,BS3: beam splitters; Mu-metal: spin rotators in both paths; <pQ and 0j: phase shifters; SR: spin rotator by a variable 
angle a; Be,: constant magnetic field. A detector counts all spin-up neutrons that leave the device via the O-beam. Another detector counts all 
neutrons that leave the device via the H-beam. 



factors. As the state vector propagates through the interferometer and the spin rotator (see Fig.fTTTi, it changes according to 

, _ / cos(a/2) zsin(a/2) \ / t* r\ ( t* r\ ( \ / e* \ 
- /sin(a/2) cos(a/2); 67 l v -r* t ) 5J { -r* t ) „\ J 6J { e<*> ) A5 

I r \ I i r \ ft -r* \ ( t -r 



r '/3.7\ r l /2,6\ J 1,5 V r ' /0,4 

-v^ 1/J) 2 , 3 ( !/^) 0>1 (' 'Oi^' "0o, 2 w ' (21) 

where the subscripts i,j refer to the pair of elements of the eight-dimensional vector on which the matrix acts. Reading back- 
wards, the first pair of matrices in Eq. d2TT i represents beam splitter BSO, the second pair the mu-metal (a spin rotation about 
the j-axis by %j\ and — 7T/4, respectively), the third and fourth pair beam splitters BS1 and BS2, respectively, the fifth pair the 
phase shifters, the sixth pair beam splitter BS3, and the last matrix represents the spin rotator SR. 

From Eq. d2"Tl ), it follows that the probability to detect a neutron with spin up in the O-beam is given by 

Po(a, X ) - |^' 3 . t | 2 = r/; 2 [l+cos(a + ^)], (22) 

where % = </>o — 0i . From Eq. (l22l it follows that the correlation Eo(a,x) is given by l53ll 

( j s Po(«,*)+Po(« + *,* + *^ (23) 
Po (a,x)+Po(a + 7c,x + n)+po(a + n,X)+Po(a,X + n) 

independent of the reflection R = \r\ 2 = 1 — T of the beam splitters (which have been assumed to be identical). Repeating the 
calculation for the probability of detecting a neutron in the H-beam shows that Eu(a,%) = 0, independent of the direction of 
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FIG. 12: Correlation E(a,x) between spin and path degree of freedom as obtained from an event-by-event simulation of a single-neutron 
interferometry experiment which shows violations of a Bell inequality 15311 . Solid surface: E(a,x) = cos(a + X) predicted by quantum 
theory; solid circles: simulation data. The lines connecting the markers are guides to the eye only. Model parameters: reflection R = 0.2 and 
J — 0.99. For each pair (a,x),N = 10000 particles were used to determine each of the four counts N(a,x)-, N(a + K,x + Tt), N(a,X + x), 
and N(a + K,x + ft) that appear in Eq. J23t . 



the spin. Note that if the mu-metal would rotate the spin about the jt-axis instead of about the y-axis, we would find Eo(a,x) = 
cos a cos x, a typical expression for a quantum system in a product state. 

The fact that Eo(a,x) = cos(a + %) implies that the state of the neutron cannot be written as a product of the state of the 
spin and the phase. In other words, in quantum parlance, the spin- and phase-degree-of-freedom are entangled 11531 15511 . In this 
context, it is customary to form the Bell-CHSH function ll56l l57ll 

S = S(ay,Zi_,ai,%i) = Eo{a x ,X\)+Eo{a\:X2)-Eo{a2,X\)+Eo{a 2 ,X2), (24) 

for some set of experimental settings oti, %\, (Xa, and Xi- If the quantum system can be described by a product state, we must 
have \S\ < 2. Therefore, if experiment shows that \S\ > 2, it is impossible to interpret the experimental result in terms of a 
quantum system in the product state |Ht]. If (X\ — Q,X\— OC2 = 7i/2, and X2 = it/4, then S = 2\[2, the maximum value 
allowed by quantum theory f58ll . 

The single-neutron interferometry experiment yields the count rate N(a,x) for the spin-rotation angle a and the difference % 
of the phase shifts of the two different paths in the interferometer 15311 . Following Ref. [53], the correlation E(a,x) is defined 
by 

£ , , = N(a,x) +N{a + K,x + n)-N{a,x + n) -N(a + n,x) 
N{a,x)+N{a + n,x + x)+N{a,x + x)+N{a + n,xy 

and, if quantum theory describes this experiment, we expect that E(a,x) ~ Eo(a,x)- Experiments show that S > 2 ll53ll54Tl . 



B. Event-by-event model: realization 

The components that constitute the interferometer have been described in Section [Til] In the following, we specify the action 
of the remaining components, namely the magnetic -prism polarizer (not shown), the mu-metal spin-turner, the spin-rotator and 
spin analyzer. 

Magnetic-prism polarizer: this device takes as input a neutron with an unknown magnetic moment and produces a neutron with 
a magnetic moment that is either parallel (spin up) or anti-parallel (spin down) with respect to the z-axis (which by definition 
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FIG. 13: Same as Fig. [T2] except that the reflection R = 0.2 and y = 0.55. The differences between the quantum theoretical results and the 
simulation data are due to the choice y = 0.55. The event-based simulation reproduces the exact results of quantum theory if y — > 1 _ (data not 
shown). 



is parallel to the guiding field B). In the experiment, only a neutron with spin up is injected into the interferometer. Therefore, 
to simplify matters a little, we assume that the neutrons that enter the interferometer all have spin up. This assumption is easily 
incorporated in the procedure that creates the initial message by simply setting 0=0 (see Eq. (fl}). 

Mu-metal spin-turner: the action of this component is to rotate the magnetic moment of the neutron by n/2 {—n/2) about the y- 
axis, depending on whether the neutron was transmitted (reflected) by BSO. The processor that accomplishes this is very simple. 
It takes as input the direction of the magnetic moment, represented by the message y and performs the rotation y «— e mc ' 4 y 
which corresponds to a rotation about the y-axis by n/2. We emphasize that we use Pauli matrices as a convenient tool to express 
rotations in 3D space, not because in quantum theory the magnetic moment of the neutron is represented by spin- 1/2 operators. 

Spin-rotator: the action of this component is to rotate the magnetic moment of a neutron by an angle a about the x-axis. It 
changes the message according to y «— e' aa ^ 2 y. 

Spin analyzer: this component selects neutrons with spin up, after which these neutrons are counted by a 100% efficient 
detector. The simplest algorithm that performs this task is to project the magnetic moment on the z-axis and send the neutron to 
the detector if the projected value exceeds a pseudo-random number $%. 

Detector: We simply count all neutrons that appear in the O- and H-beam. 



C. Simulation results 



In Fig. Q~2] we present results for the correlation Eq. d25l l as obtained from event-by-event simulations of the experimental 
setup depicted in Fig.QT] assuming that the experimental conditions are very close to ideal. For the ideal experiment, quantum 
theory predicts that E(a, x) = cos(a + % ) (represented by the solid surface in Fig. [121) and as shown by the markers in Fig. [12] 
disregarding the small statistical fluctuations, there is close-to-prefect agreement between the event-by-event simulation data and 
quantum theory. 

The real experiment suffers from unavoidable imperfections, leading to a reduction and distortion of the interference pat- 
terns [53]. In an event-by-event approach, it is easy to incorporate mechanisms for different sources of imperfections by mod- 
ifying or adding rules. After all, we can manipulate each individual event. However, to reproduce the available data, this is 
not necessary because, as before, we can use the parameter y to control the deviation from the quantum theoretical result. In 
particular, we can use the parameter y to fit the simulation results to the experimental data for the value of S. 
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For instance, taking R = 0.2 and 7 = 0.55, the simulation (see Fig. \13[ yields Smax = 2.05, in excellent agreement with the 
value 2.052 ± 0.019 obtained in experiment [53]. For R = 0.2 and 7 = 0.67, the simulation yields S mux = 2.30, in very good 
agreement with the value 2.291 ± 0.008 obtained in a similar, more recent experiment lf54Tl . 

D. Discussion 

An Einstein-Podolsky-Rosen-Bohm (EPRB) experiment can be used to test for violations of a Boole-Bell-type inequality but 
not all experiments that test for violations of a Boole-Bell-type inequality are EPRB experiments. Essential features of an EPRB 
(thought) experiment are that (1) a source emits pairs of particles with properties of which at least one is correlated, (2) as the 
particles leave the source, they no longer interact (but the correlations of their properties do not change), (3) the properties of the 
particles are determined by two spatially separated analyzers which do not communicate with each other and have settings that 
may be changed independently and randomly, (4) all particles leaving the source are analyzed and contribute to the averages and 
correlations. Clearly, the neutron interferometry experiment which we have discussed in this section is not an EPRB experiment. 
It only satisfies the fourth criterion if we disregard the neutrons that are transmitted by BS1 or BS2. 

The violation of the Bell-CHSH inequality observed in the neutron interferometry experiment demonstrates that it is possible 
to create a correlation between the path and the magnetic moment of the neutron, although there is no direct "interaction" 
between the two. If we interpret the outcome of this experiment in terms of quantum theory, the observed violation of the Bell- 
CHSH inequality implies that it is impossible to describe the outcome of the experiment in terms of a product state of path and 
spin states. Hence, the system must be described by an entangled state. On the other hand, a classical, Einstein-local and causal, 
event-by-event process can also reproduce all the features of the entangled state. Hence, not too much significance should be 
attached to the latter. 

For completeness, we mention that if we pick the angle % randomly from the same finite set of predetermined values used 
to produce Fig. [T2l an event-based simulation with 7 = 0.99 yields (within the usual statistical fluctuations) the correlation 
E(a,x) ~ (1 /2)cos(a +X)> which does not lead to a violation of a Bell-type inequality (data not shown). Thus, if the neutron 
interferometry experiment could be repeated with random choices for the position of the phase shifter (x), and the experimental 
results would show a significant violation of a Bell-type inequality, the event-based model that we have presented here would be 
ruled out. 



VI. COHERENT MANIPULATION OF THE NEUTRON SPIN, PHASE AND ENERGY 

In the language of quantum theory, the neutron interferometry experiment described in Ref. ll59ll demonstrates that using 
magnetic resonance techniques, it is possible to create entanglement between the spin, path and energy degree of freedom of the 
neutron. 

The set-up of the single-neutron interferometry experiment for the observation of interference between the phases induced 
by a phase shifter (PS) and by two radio-frequency (RF) fields is shown in Fig. [14] (see also Fig. 1 in Ref. f59ID . Neutrons 
with their spin polarized along the static magnetic field B = Be z impinge on the interferometer. The radio-frequency coil (RF1) 
within one path of the skew-symmetric neutron interferometer is tuned such that it flips the spin from up to down and induces 
a phase shift (<f>i). Neutrons that leave the interferometer in the H-beam are counted by a neutron detector. Neutrons that leave 
the interferometer in the O-beam pass an accelerator coil (AC) that compensates for the differences in times of flight of the two 
different paths in the interferometer. These neutrons pass through a second radio-frequency coil (RF2) which is placed in a 
region where the static magnetic field is B/2 and which is driven by the half frequency of that of RF1. RF2 flips the spin and 
induces a phase shift (^2)- Subsequently, neutrons with their spin up are counted by the combination of a n/2 spin rotator (SR), 
an analyzer and a neutron detector. 

A. Quantum theory 

In Fig. [15] we show the diagram that corresponds to the experiment ll59ll . Spin-polarized neutrons impinge on beam splitter 
BS0. A constant magnetic field (Be z ) defines the spin-up direction. The Hamiltonian that describes the spin of the neutron as it 
moves through the radio-frequency coil RF1 reads 

Hi = -Be-S z -Bie y S y cos(cot + (pi), (26) 

where B and Bi are the static and radio-frequency (RF) fields and CO and cpi denote the angular frequency and phase of the RF 
field, respectively. The wave function of the neutron spin, denoted by |*P(f )) satisfies the time-dependent Schrodinger equation 
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FIG. 14: A single-neutron interferometry experiment for the observation of interference between the phases induced by a phase shifter (PS) 
and by two RF fields [59]. BS0,...,BS3: beam splitters; MFP: magnetic field prism; B = fie-: constant magnetic field; PS: phase shifter ; RF1, 
RF2: radio-frequency spin-flippers; AC: accelerator coils; ST: spin turner. 



(TDSE) 



(27) 



where from now on, we use units such that ft = 1. Writing = e'^'+Pi^V^ 1 ' 5 "*/ 2 !*!^?)), and imposing the resonance 

condition B = co, the TDSE for |4>(f )) reads id \®(t)) jdt = H[ (t)\®(t)) where 



H[(t) = - 1 - 



B t B t 

5 v cos2(<af + tpi) - S x sm2(cot + q>\ ) cos — — 5 z sin2(©? + <p{) sin — 



(28) 



showing that H[(t) oscillates with frequencies 2(0 and 2co ±B\ 
can be neglected, as is usually assumed in NMR/ESR theory |6C 
its time evolution is given by the unitary matrix 



2. We assume that the effect of these oscillations on the spin 
1621 1 . As the neutron spin passes through the spin flipper RF1, 



Ui(t) 



(29) 



where t\ is the time during which the neutron experiences the RF field of RF1. Adjusting B\ such that B\t\/2 = %, the RF field 
changes the neutron spin from up to down (see Fig. [15]) and changes the phase by (Ot\ + <p\. 

As the neutron leaves BS3, its spin can be up or down. All neutrons that leave BS3 via the H-beam are sent to a detector. 
Neutrons with spin up (down) that are transmitted (reflected) by BS3 fly through the spin-flipper RF2, operating at a>2 — co/2 at 
resonance with the static field B/2 that is present in region 2. The two RF-spin-flippers and the static fields in regions 1 and 2 act 
as an interferometer for the energy of the neutron [59]. Finally, the neutrons pass a spin rotator SR which rotates the magnetic 
moment of the neutron about the x-axis by 7z/2, mixing the spin-up and spin-down components. 
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FIG. 15: Quantum theoretical model of the single-neutron interferometry experiment that demonstrates the entanglement of the spin, path and 
energy degrees of freedom 1 59] . Polarizer spin-up neutrons are injected in an interferometer. BS0,...,BS3: beam splitters; tj>o and phase 
shifters; RF1, RF2 spin flippers: radio frequency coils, tuned to the resonance frequencies of the neutron spin in the static magnetic fields B 
and B/2, respectively; SR: spin rotator (magnet causing a magnetic moment to rotate about the x-axis by Tt/2). A detector counts the spin-up 
neutrons in the O-beam. Another detector counts all neutrons in the H-beam. 



As the state vector propagates through the interferometer and the spin analyzer, see Fig. |T5j it changes according to 



1/V2 i/y/2 \ ( c«2+ft)/2\ / t * r\ ft* 

i/y/l l/V2j 6J { -^1^)12 o ) 78 [ -r* t) 5 \ -r* t / ^ 

t* r\ ( t* r\ ft -r* \ ( t -r* \ ( 1 



r 



-r* t I \ —r* t I \ r t* I \ r t* I \ e imt i , 

J 3,7 V / 2.6 V / 1.5 V / 0,4 V U e / 4,5 

e i(m+9i)/2 \ / e i<t>i o \ /e* \ ft -r*\ ft -r* 

- e «'(- 1+ ^)/2 o ) \ ) 2i3 [ e^) 0t \r f ) [r t* , 



m- oo) 



Reading backwards, the matrices in Eq. (130b represent beam splitter BSO (first pair of matrices), the phase shifters inducing a 
relative phase shift % = <po — 0i (second pair) , the RF-spin-fiipper RF1 (fifth matrix), a phase shift (sixth matrix) applied to 
the spin-down neutron (f| is the time during which a neutron appearing in the O-beam has spin down) beam splitters BS1, BS2 
and BS3 (matrices 7 through 12), the RF-spin-flipper RF2 (matrix 13), and finally the spin rotator SR (a spin rotation about the 
x-axis by 7C/2). The analyzer (not included in Eq. d30l l) sends the spin-up neutrons in the O-beam to the detector. 
From Eq. ( l30l l. it follows that the probability to detect a neutron with spin up in the O-beam is given by 

Po((p,X) = \ x i" 3 /-^TR 2 {l+sm[x + (Ot l + (t 2 -t l )(o/2 + ( P ]} 1 (31) 

where <p = ^ — (pi /2, showing that the spin (<pi,(P2), path (x), and energy (flWp degree-of-freedom are entangled and can be 
manipulated independently. In the experiment 15911 . the phases induced by the guiding fields in regions 1 and 2 and the phase 
+ ( f 2 — t\)(o/2 were compensated for by a tunable accelerator coil with a static magnetic field along the z-direction lf59tl . 
Putting (Of), + (t 2 -ti)C0/2 = in Eq. ([31} yields 

Po ((p,X) = \%, t \ 2 = TR 2 {l+sm{x + (p}}. (32) 
The probability to find a neutron leaving the device via the H-beam reads 

PB = |^, t | 2 + |^, ; | 2 = J R(r 2 + J R 2 ), (33) 
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FIG. 16: Normalized counts N(<p,x) as obtained from an event-by-event simulation of a single-neutron interferometry experiment that employs 
two radio-frequency fields to manipulate the energy of a single neutron [59]. Solid surface: probability po((p,%) = (\ — R)R 2 [l + sin(^; + <p)] 
to observe neutrons in the O-beam, as predicted by quantum theory; solid circles and open squares: simulation data for the normalized neutron 
counts in the O- and H-beam, respectively. The lines connecting the markers are guides to the eye only. Quantum theory predicts the probability 
to observe a neutron in the H-beam to be R(R 2 + (1 — R) 2 ) = 1 /4. Model parameters: reflection R = 0.5, y = 0.99; number of particles per 
data point N = 10000. 



which obviously does not depend on %, <pi, t\, q>2, or t2- 



B. Event-by-event model: realization 

The components that constitute the interferometer have been described in Section|III]and the components that manipulate the 
neutron spin (in the rotating frame in the case of the RF flippers) have been described in Section IVBl Therefore, we simply use 
these components without modification. Also the simulation procedure is the same as before. The source sends a particle to BS0 
and it is not until that particle has left the interferometer or has been detected in either the O- or H-beam that the source will send 
the next particle. 



C. Simulation results 



In Fig. [16] we present results for the normalized particle count N((p,x) in the O-beam (solid circles) and H-beam (open 
squares) in the case of the ideal experiment. The normalization consists of dividing the actual count by the number of particles 
N generated by the source. The quantum theoretical expression for the normalized O-beam count is given by Eq. ( l32l l and is 
represented by the solid surface in Fig. [16] According to Eq. (l33l . quantum theory predicts that the normalized H-beam count is 
independent of % and <p. From Fig. [16] it is obvious that the event-based simulations reproduces all the features of the quantum 
theoretical description of this experiment. 

The effect of changing the reflection from R = 0.5 (Fig.[T6ll to R = 0.2 in combination with reducing the parameter y is shown 
in Fig.[l7] As expected from Eqs. (f32t and d33l . the normalized counts are reduced and, because y = 0.5 instead of being close 
to one, the simulation data deviate (slightly) from the quantum theoretical results (solid surface) Eq. (l32l . 

Our conclusion is that the event-based particle-only model reproduces the results of quantum theory for a neutron interferom- 
etry experiment which, in the language of quantum theory, exhibits entanglement ll59ll . 
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FIG. 17: Same as Fig. Q7d] except that R = 0.2 and y = 0.5. Quantum theory predicts the probability to observe a neutron in the H-beam to 
be R(R 2 + ( 1 — R) 2 ) ~ 0. 19. Differences between the quantum theoretical results and the simulation data are due to the choice y = 0.5. The 
event-based simulation reproduces the exact results of quantum theory if y — > 1~ (data not shown). 



An important question is whether the event-based approach leads to new predictions that may be tested experimentally and 
may therefore reveal new physics. As demonstrated in this paper, and in our earlier work on (quantum) optics experiments (for 
a review see Ref. ifTTll . in the stationary state (meaning after processing many events) the event-based simulation reproduces the 
statistical distributions of quantum theory. 

Therefore, new predictions can only appear when the event-based simulation is operating in the transient regime, before the 
processors reach their stationary state. Optics experiments with a Mach-Zehnder interferometer and two-beam interference that 
may be able to address this question have been discussed in Ref. 112 ill and Ref. 11911 . respectively. Neutron interferometry is 
well-suited to address this issue because the neutrons can be detected with almost 100% efficiency and a relatively low flux of 
neutrons facilitates labeling each detection event by a time stamp. Therefore, it is feasible to collect detailed information about 
each neutron, which can then be analyzed further. 



In these experiments, interference patterns were recorded such that in a fixed time frame, the sum of neutron counts in the O- 
beam over all chosen settings of the phase shifts is approximately constant (and relatively small) |2j,|63|,|64|]. This is accomplished 
by increasing the rotation angle % of the phase shifter at regular time intervals. 

Adopting the same procedure as in experiment, the event-based simulation yields the results presented in Fig. [18] The 
simulation data show the same qualitative features as found in experiment lE l63ll . 

It is instructive to inquire why the event-based processor is able to reproduce all these features, even though it operates far 
from a stationary-state regime. The main reason for this can be traced back to the speed with which the DLM responds to a 
change in the input messages (see section|II|. In Section uTCl we explained that y controls the speed and accuracy with which 
the event-based processor responds to a change in the ratio of the input events on its input ports. In the case of the neutron 
interferometer (see Fig. [TJ, the ratio of the number of neutrons that travels from BS 1 to BS3 and from BS2 to BS3 is independent 
of the phase shift %. Therefore, as far as the response time to a change is concerned, only the message content matters. However, 
the event-based processor described in Section HTCl is constructed such that the last messages that were delivered at input ports 
and 1 are stored in the DLM. Hence, the response to a change in one of the messages is immediate and without introducing 
errors. 



VII. TIME-DEPENDENT PHENOMENA 



A. Low-counting-rate experiments 
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FIG. 18: Phase shift dependence extracted from interference patterns taken with a small number of neutrons, see also Fig. 5 in Ref.l63l The 
number of the detected neutrons in the O-beam, summed over the five different settings of the phase shift and averaged over 30 independent 
simulations is (N). Solid squares: Average counts of 30 simulations; solid line: least square fit of the average counts to a sinusoidal func- 
tion; solid circles: counts obtained from individual runs with approximately (N) detected neutrons; error bars: one standard deviation. All 
simulations where carried out with R = 0.2 and y = 0.5. 



B. Time-dependent beam blocking 

Summhammer's neutron interferometry experiment [65] with a shutter which, upon detection of a neutron, randomly blocks 
one of the paths through the interferometer provides what is perhaps the most stringent test of the event-based model which we 
present in this paper. 

The experimental setup is sketched in Fig. [19] Depending on the state of the Cd metal shutter, neutrons transmitted by BS0 
are blocked. Neutrons transmitted by BS1 or BS2 leave the interferometer and do not contribute to the neutron counts in the O- 
or H-beam. Upon detection of a neutron, the shutter may change its state with probability 1/2. Detection events are labeled by 
the state of the shutter. 

The experimental data show the following features (see Fig. 3 of Ref. i65l0 : 

1. If the shutter is kept closed, the relative frequency of the O-beam detection events does not depend on the phase shift %, 
that is there is no interference pattern. 

2. If the shutter is kept open, the relative frequency of the O-beam detection events shows the sinusoidal dependence on %, 
the characteristic signature of interference. 

3. If the state of the shutter is allowed to change according to the procedure described earlier, the relative frequency of the 
O-beam detection events conditioned on the state of the shutter exhibits a sinusoidal dependence on % when the shutter 
was open and no dependence on % when the shutter was closed. 

4. The maximum of the relative frequency of the O-beam detection events when the shutter was open is approximately equal 
to the relative frequency of the O-beam detection events when the shutter was closed. This maximum is approximately 
0.43. 

5. The visibility of the relative frequency of the O-beam detection events when the shutter was open is about 0.4. 
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FIG. 19: Layout of Summhammer's neutron interferometry experiment with a shutter 1 65] . BS0,...,BS3: beam splitters. Neutrons which are 
transmitted by BSO may be blocked by a piece of Cd metal, depending on the state of the shutter. Phase shifter: aluminum foil. Neutrons that 
are transmitted by BS1 or BS2 leave the interferometer and do not contribute to the neutron counts in the O- or H-beam. For each detected 
neutron, the state of the shutter changes with probability 1/2 1 65]. The detection events are labeled by the state of the shutter. 




FIG. 20: Event-based simulation results of Summhammer's neutron interferometry experiment with a time-dependent shutter 165], see Fig.l 191 
The relative frequency is the neutron count in the O-beam divided by the sum of the counts in the O- and H-beam. Open (solid) circles: relative 
frequency of events recorded with the shutter open (closed), a) Shutter closed; b) Shutter open; c) For each neutron detected in the O- or 
H-beam, the state of the shutter changes with probability 1/2. Simulation parameters: number of incident neutrons N = 250000, reflection 
R — 0.4, and y = 0.12. In case (c), the total number of events (per value of x) registered in the O- and H-beam is approximately 60000, as in 
experiment 1 65] . 



Originally conceived to test a nonergodic interpretation of quantum theory l66l - l68ll . the experimental results were interpreted 
as being in accordance with the Copenhagen interpretation and ruling against predictions based on the nonergodic interpretation 
of quantum theory l65l l69ll . Citing Summhammer about the latter, "This prediction expects that in a neutron interferometer 
successive neutrons interact with each other through a hypothesized memory to which each neutron contributes a little, such that 
the quantum mechanical interference phenomenon only arises when sufficiently many neutrons have passed the interferometer 
in a constant experimental condition". 

As explained in Section HU the key feature of the event-based model of a beam splitter is that it can adapt to changes of the 
input data, in other words, it can learn. Obviously, learning requires some form of memory, which in our case consists of the x 
and Y registers, see Fig. [2] As Summhammer's experiment rules against the nonergodic interpretation of quantum theory B65l 
and memory is a key feature of this interpretation, one might expect that this experiment rules out the event-based approach as 
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FIG. 21: Same as Fig.[20]using the simplest approximation to account for the incoherence of the incident neutron beam (see text). Simulation 
parameters: number of incident neutrons N = 500000, reflections R\ = 0.2, R2 = 0.9, weights W\ =20/21, W2 = 1/21, and 7 = 0.12. In case 
(c), the total number of events (per value of x) registered in the O- and H-beam is approximately 60000, as in experiment 16511 . 



well. However, as we now show, this is not the case. 

For all experiments considered in this paper, it is essential to include into the event-based model all the essential aspects of the 
real experiments, not just those that are considered relevant on the level of idealized thought-experiments. In Summhammer's 
experiment, the motion of the shutter induces vibrations in the silicon crystal ll65l l69ll . In the absence of concrete knowledge 
about this effect, it seems very difficult to model in detail how the opening and closing of the shutter affects the vibrating crystal. 
Therefore, we adopt a pragmatic approach to mimic the effect of the shutter motion on the interferometer. Of course, we could try 
different update rules for the event-based processors but even that is not necessary. From computer experiments, we found that 
it is sufficient to reset the internal vectors x to zero each time the shutter closes. Alternatively, we can set their respective values 
of 7 to zero each time the shutter closes and reset 7 of a particular processor to its specified value after it has sent out a message. 
Both these modifications have a simple physical interpretation in terms of shaking the crystal and temporarily destroying the 
coherence in (parts of) the silicon crystal. Obviously, the ability to easily incorporate the effect of such processes is a powerful 
feature of a discrete-event simulation approach. 

A first set of simulation results is shown in Fig. [20] It is clear that the event-based simulation reproduces the main features 
(see items 1-3 above) of Summhammer's experimental results [65]. What is still missing is the quantitative agreement with 
Summhammer's data. In fact, we have found it to be impossible to use the parameters R and 7 to achieve good numerical agree- 
ment. This simply means that we should consider moving away from the effective, average reflection coefficient characterization 
of the beam splitter, in agreement with the theory given in Ref. |2]. 

Instead of using one reflection coefficient R, let us try to use only two reflection coefficients R\ and R2 (for BS0,...,BS3). Two 
corresponding weights W\ and W2 determine the frequency with which the individual neutron "sees" one of the two reflection 
coefficients. This approximation can be viewed as a minimalistic substitute of the incoherent averaging over the Pendellosing 
oscillations and misset angle in the wave theoretical treatment of the perfect crystal neutron interferometer (see Chapter 10 of 
Ref. (2[]). As Fig.|2T|demonstrates, this minimal extension suffices to reach quantitative agreement with Summhammer's data, 
that is the simulation reproduces all five features listed above. 



Quantum theoretical treatment 



It is instructive to scrutinze the statement that the experimental results are in accordance with the Copenhagen interpretation 
of quantum theory l65ll . Adopting the effective, one-parameter model of the beam splitter, in the case that the shutter is open, 
quantum theory predicts that the probabilities to observe a neutron are given by (see Eq. ( fToT l) 

p open = 2TR 2 (1+vcosx) 

p° pen = R(T 2 +R 2 -2vRTcos X ) 

open _ T 2 

Pbs\ — 1 

Pbs? = TR i 04) 

where and p^ 2 " are tne probabilities that the neutrons are transmitted by BS1 and BS2, respectively. Anticipating for 
the observation that the visibility of the interference fringes is (much) less than one, we have introduced the visibility v as an 
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adjustable parameter. On the other hand, if the shutter is closed we have 

dosed t 1 d2 
Po —IK 



dosed n3 

Ph k 

.dosed 

PCd — 1 

pt°s? d = RT, (35) 

where /?cd Sed represents the probability that the neutron is absorbed by the Cd shutter. 

From Eqs. (f34b and d35l ). it follows that the relative frequencies to observe neutrons in the O-beam are given by 

open 

/° Pen (Z) - o P en° open =2TR{\+ VCOSj) , (36) 

Po +Ph 

„closed 

^•closed £ O y (37) 

J „closed i ...dosed ' ^ ' 

Po + Ph 

The experiment shows that max z / open (^) = y closed (see (4) in the list above) hence 2(1 + v)TR = T or, using the observation 
that v ss 0.4 (see (5) in the list above), R = 1/2(1 + v) ss 0.36 such that / closed « 0.64, which is incompatable with the experi- 
mental observation that y closed « 0.43 (see (4) in the list above). Therefore, the experimental data rules out the two-parameter 
(R, v) quantum model of this interferometry experiment. However, just as we did in the case of the event-based model for this 
experiment, we may generalize the model summing over incoherent processes. As before, for simplicity, we consider a model 
with only two reflections R\ and R2 occurring with probabilities W\ and W2 = 1 — W\, respectively. 
Instead of Eq. d37l ). we now have 

r ^ = 2 w.t iR i + w 2 t 2 rI 

J WiRi+W 2 R 2 
losed = W l T l R 2 l +W 2 T 2 Rl 
W x R\+W 2 Rl ' 

and imposing the conditions max^/ open (j) = y closed = g yields 

2(l+v)(W l R] + W 2 R2) = W1R1+W2R2, 

WiT l R 2 l +W 2 T 2 Rl = g{WiR} + W 2 Rl). (39) 

These equations can be solved in closed form, the non-trivial solution reading 



^2 



l-2(l+v)/?i' 



Wl = (l- g -^)(-2v + 2 g (v + l)-l) 

8(y+l) 3 /?{-12(v+l) 2 ^+6(v+l)^-2(g + ( g -l)v)/? 1 -( g -l)(-2v + 2 g (y+l)-l)' 

Note that the solutions are bounded by < R\ < 1, < R 2 < 1, < v < 1, < g < 1 and < W\ < 1. 

Taking for instance R\ = 0.2, v = 0.4 and g — 0.43, we obtain R2 = 0.84 and W\ = 1 — W2 = 0.93, which are all reasonable 
numbers, the values of R\,R2,Wi and W2 being rather close to the numbers that were used in the event-based simulation (see 
Fig- ED- In conclusion, it is clear that the four-parameter (Ri,R 2 ,Wi,v) quantum model is compatible with the experimental 
data. Note that this compatibility is not due to some unique feature of quantum theory but merely results from adding, with 
appropriate weights, the results of two independent experiments performed under different conditions. 

In spite of the fact that quantum theory can be used to describe the outcome of Summhammer's experiment, the "mystery" 
alluded to in the introduction remains. First, it is mysterious how the experimental apparatus can "know" the expresssions of the 
probabilities Eqs. (T34t and ( 1331 ) before the very first neutron has passed through it and a decision about the state of the shutter 
was taken. Second, quantum theory postulates that a detection event corresponds to a certain value (one out of four in this case) 
of a random variable with a probability distribution given by either Eq. d34i l or Eq. d35l l but (like probability theory in general) 
is silent about the process by which these random values are actually realized. There are no such mysteries in the event-based 
approach as it provides a complete prescription of how individual (detection) events are to be generated. 
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VIII. SUMMARY AND OUTLOOK 

In this paper, we have demonstrated that the event-based approach, originally introduced in Ref. J^, [l(l [l5ll to simulate 
quantum optics experiments, can also be applied to neutron interferometry experiments. Our approach gives a detailed, mystery- 
free, particle-only description of interference and entanglement, as observed in neutron interferometry experiments J2[] and does 
not suffer from the quantum measurement problem, simply because the discrete events, such as the detection of a neutron, are 
taken as the basic entities of the description. The statistical distributions which are observed in real experiments, usually thought 
to be of quantum mechanical origin, emerge from a time series of discrete events generated by causal, classical, adaptive systems. 

Conceptually, our approach may seem to have similarities to cellular autonoma modeling ll70lF7Tll . the work of 't Hooft ll72ll73ll 
or to, for instance, lattice Boltzmann modeling of fluid dynamics B74I1 which all explore the idea that simple rules, which not 
necessarily derive from classical Hamiltonian dynamics, define discrete-event processes which may lead to the (complicated) 
behavior that is observed in experiments. However, the reasoning that lead us to our simulation model is very different. 

Starting from the point of view that empirical knowledge, and the concepts created on the basis of this knowledge, derives 
from the elementary events which are registered by our senses, we explore the consequences of assuming that current scientific 
knowledge is built on the notion of discrete events and the relations between them. This is a departure from the prevailing mode 
of thinking in theoretical physics, which assumes that the definite results which we observe are signatures of an underlying 
objective reality that is mathematical in nature. While the hypothesis that such a reality exists cannot be refuted on logical 
grounds, their is no experimental evidence that supports it. Taking an indifferent stance on this issue, the urge to "deduce" the 
existence of definite results (discrete events) from a set of axioms disappears, opening a route to a mode of thinking that is 
much less constrained. Apparently, this change of paradigm facilitates the construction of simulation models which reproduce 
the experimental and quantum theoretical results of many real experiments, in particular those in which the data is recorded 
event-by-event. 

The discrete-event model that has been described in this paper can also be used to simulate neutron spin-echo experiments dH 
l75ll and the recent experiments on the uncertainty in neutron spin measurements [76]. However, neutron (and optics) experiments 
that involve diffraction/scattering cannot yet be simulated with the present model. Incorporating this feature is left for future 
research. 

We hope that our work will stimulate the design of new single-neutron experiments to explore the applicability of event-based 
modeling to physical phenomena. Specifically, to test the validity of our discrete-event modeling, it would be worthwhile to 
repeat Summhammer's experiment with the shutter under conditions that show much larger visibility. 

IX. ACKNOWLEDGEMENT 

We are grateful to Profs. H. Rauch and J. Summhammer for very stimulating discussions and Dr. H. Lemmel for making 
experimental data available to us. We are indebted to K. De Raedt for his help in solving the time-dependent beam blocking 
problem and to Profs. M. Novotny and D. Stavenga for a critical reading of the manuscript. This work is partially supported by 
NCF, the Netherlands. 



[1] A. Tonomura, The Quantum World Unveiled by Electron Waves (World Scientific, Singapore, 1998). 

[2] H. Rauch and S. A. Werner, Neutron Interferometry: Lessons in Experimental Quantum Mechanics (Clarendon, London, 2000). 

[3] J. C. Garrison and R. Y. Chiao, Quantum Optics (Oxford University Press, Oxford, 2009). 

[4] D. Home, Conceptual Foundations of Quantum Physics (Plenum Press, New York, 1997). 

[5] L. E. Ballentine, Quantum Mechanics: A Modern Development (World Scientific, Singapore, 2003). 

[6] A. Leggett, in The Lessons of Quantum Theory: Niels Bohr Centenary Symposium, edited by J. de Boer, E. Dal, and O. Ulfbeck (Elsevier, 

Amsterdam, 1987), pp. 35 - 58. 
[7] A. Allahverdyan, R. Balian, and T. Nieuwenhuizen, arXiv: 1 107.2138 (201 1). 

[8] R. P. Feynman, R. B. Leighton, and M. Sands, The Feynman Lectures on Physics, Vol. 3 (Addison- Wesley, Reading MA, 1965). 

[9] K. De Raedt, H. De Raedt, and K. Michielsen, Comp. Phys. Comm. 171, 19 (2005). 
[10] H. De Raedt, K. De Raedt, and K. Michielsen, J. Phys. Soc. Jpn. Suppl. 76, 16 (2005). 
[11] K. Michielsen, F. Jin, and H. De Raedt, J. Comp. Theor. Nanosci. 8, 1052 (201 1). 
[12] H. Rauch, W. Treimer, and U. Bonse, Phys. Lett. A 47, 369 (1974). 
[13] Y. Hasegawa and H. Rauch, New J. Phys. 13, 115010 (2011). 

[14] G. Kroupa, G. Bruckner, O. Bolik, M. Zawisky, M. Hainbuchner, G. Badurek, R. Buchelt, A. Schricker, and H. Rauch, Nucl. Instrum. 

Methods Phys. Res. A. 440, 604 (2000). 
[15] H. De Raedt, K. De Raedt, and K. Michielsen, Europhys. Lett. 69, 861 (2005). 
[16] S. Zhao, S. Yuan, H. De Raedt, and K. Michielsen, Europhys. Lett. 82, 40004 (2008). 



30 



[17] K. Michielsen, S. Yuan, S. Zhao, F. Jin, and H. De Raedt, Physica E 42, 348 (2010). 
[18] F. Jin, S. Zhao, S. Yuan, H. De Raedt, and K. Michielsen, J. Comp. Theor. Nanosci. 7, 1771 (2010). 
[19] F. Jin, S. Yuan, H. De Raedt, K. Michielsen, and S. Miyashita, J. Phys. Soc. Jpn. 79, 074401 (2010). 
[20] F. Jin, H. De Raedt, and K. Michielsen, Commun. Comput. Phys. 7, 813 (2010). 

[21] K. Michielsen, Th. Lippert, M. Richter, B. Barbara, S. Miyashita, and H. De Raedt, J. Phys. Soc. Jpn. 81, 034001 (2012). 
[22] P. Grangier, G. Roger, and A. Aspect, Europhys. Lett. 1, 173 (1986). 

[23] V. Jacques, E. Wu, F. Grosshans, F. Treussart, P. Grangier, A. Aspect, and J.-F. Roch, Science 315, 966 (2007). 
[24] P. D. D. Schwindt, P. G. Kwiat, and B.-G. Englert, Phys. Rev. A 60, 4285 (1999). 

[25] V. Jacques, E. Wu, T. Toury, F. Treussart, A. Aspect, P. Grangier, and J.-F. Roch, Eur. Phys. J. D 35, 561 (2005). 
[26] S. Zhao and H. De Raedt, J. Comp. Theor. Nanosci. 5, 490 (2008). 

[27] I. N. Agafonov, M. V. Chekhova, T. S. Iskhakov, and A. N. Penin, Phys. Rev. A 77, 053801 (2008). 
[28] K. Michielsen, K. De Raedt, and H. De Raedt, J. Comput. Theor. Nanosci. 2, 227 (2005). 
[29] A. Aspect, P. Grangier, and G. Roger, Phys. Rev. Lett. 49, 91 (1982). 
[30] A. Aspect, J. Dalibard, and G. Roger, Phys. Rev. Lett. 49, 1804 (1982). 

[31] G. Weihs, T. Jennewein, C. Simon, H. Weinfurther, and A. Zeilinger, Phys. Rev. Lett. 81, 5039 (1998). 
[32] K. De Raedt, K. Keimpema, H. De Raedt, K. Michielsen, and S. Miyashita, Euro. Phys. J. B 53, 139 (2006). 
[33] H. De Raedt, K. De Raedt, K. Michielsen, K. Keimpema, and S. Miyashita, J. Phys. Soc. Jpn. 76, 104005 (2007). 
[34] K. De Raedt, H. De Raedt, and K. Michielsen, Comp. Phys. Comm. 176, 642 (2007). 

[35] H. De Raedt, K. De Raedt, K. Michielsen, K. Keimpema, and S. Miyashita, J. Comp. Theor. Nanosci. 4, 957 (2007). 

[36] H. De Raedt, K. Michielsen, S. Miyashita, and K. Keimpema, Euro. Phys. J. B 58, 55 (2007). 

[37] S. Zhao, H. De Raedt, and K. Michielsen, Found, of Phys. 38, 322 (2008). 

[38] B. Trieu, K. Michielsen, and H. De Raedt, Comp. Phys. Comm. 182, 726 (2011). 

[39] A. W. Overhauser and R. Colella, Phys. Rev. Lett. 33, 1237 (1974). 

[40] R. Colella, A. W. Overhauser, and S. A. Werner, Phys. Rev. Lett. 34, 1472 (1975). 

[41] T. Jenke, P. Geltenbort, H. Lemmel, and H. Abele, Nat. Phys. 7, 468 (201 1). 

[42] T. Unnerstall, Phys. Lett. A 151, 263 (1990). 

[43] M. Bom and E. Wolf, Principles of Optics (Pergamon, Oxford, 1964). 

[44] A. Taflove and S. Hagness, Computational Electrodynamics: The Finite-Difference Time-Domain Method (Artech House, Boston, 2005). 

[45] H. Rauch and M. Suda, Phys. Stat. Sol. A 25, 495 (1974). 

[46] M. A. Home, Physica B+C 137, 260 (1986). 

[47] H. Lemmel and A. G. Wagh, Phys. Rev. A 82, 033626 (2010). 

[48] H. Rauch and J. Summhammer, Phys. Lett. A 104, 44 (1984). 

[49] J. Summhammer, H. Rauch, and D. Tuppinger, Phys. Rev. A 36, 4447 (1987). 

[50] H. Rauch and J. Summhammer, Phys. Rev. A 46, 7284 (1992). 

[51] P. Kaloyerou and H. Brown, Physica B: Condensed Matter 176, 78 (1992). 

[52] G. R. Grimmet and D. R. Stirzaker, Probability and Random Processes (Clarendon Press, Oxford, 1995). 
[53] Y. Hasegawa, R. Loidl, G. Badurek, M. Baron, and H. Rauch, Nature 425, 45 (2003). 

[54] H. Bartosik, J. Klepp, C. Schmitzer, S. Sponar, A. Cabello, H. Rauch, and Y. Hasegawa, Phys. Rev. Lett. 103, 040403 (2009). 
[55] S. Basu, S. Bandyopadhyay, G. Kar, and D. Home, Phys. Lett. A 279, 281 (2001). 

[56] J. S. Bell, Speakable and Unspeakable in Quantum Mechanics (Cambridge University Press, Cambridge, 1993). 
[57] J. F. Clauser, M. A. Home, A. Shimony, and R. A. Holt, Phys. Rev. Lett. 23, 880 (1969). 
[58] B. S. Cirel'son, Lett. Math. Phys. 4, 93 (1980). 

[59] S. Sponar, J. Klepp, R. Loidl, S. Filipp, G. Badurek, Y. Hasegawa, and H. Rauch, Phys. Rev. A 78, 061604 (2008). 
[60] A. Abragam, The Principles of Nuclear Magnetism (Clarendon Press, Oxford, 1961). 

[61] R. Ernst, G. Bodenhausen, and A. Wokaun, Principles of Nuclear Magnetic Resonance in One and Two Dimensions, International series 

of monographs on chemistry (Clarendon Press, 1990). 
[62] J. Weil, J. Bolton, and J. Wertz, Electron paramagnetic resonance: elementary theory and practical applications (Wiley, 1994). 
[63] H. Rauch, J. Summhammer, M. Zawisky, and E. Jericha, Phys. Rev. A 42, 3726 (1990). 
[64] M. Zawisky, Master's thesis, Technische Universitat Wien (1993). 
[65] J. Summhammer, Nuovo Cimento. B 103, 265 (1989). 
[66] V. Buonomano, Nuovo Cimento B 57, 146 (1980). 
[67] V. Buonomano, Nuovo Cimento B 43, 69 (1985). 
[68] V. Buonomano and F. Bartmann, Nuovo Cimento B 95, 99 (1986). 
[69] V. Buonomano, Found. Phys. Lett. 2, 565 (1989). 
[70] E. Fredkin, Physica D 45, 254 (1990). 

[71] S. Wolfram, A New Kind of Science (Wolfram Media Inc., 2002). 

[72] G. 't Hooft, in Beyond the Quantum, edited by T. M. Nieuwenhuizen, B. Mehmani, V. Spicka, M. Aghdami, and A. Y. Khrennikov (World 

Scientific, 2007), pp. 3 - 19. 
[73] G. 't Hooft (2012), URL larXiv : 1204 . 4926vl1 

[74] S. Succi, The Lattice Boltzmann Equation for Fluid Dynamics and Beyond (Clarendon Press, Oxford, 2001). 

[75] S. V. Grigoriev, W. H. Kraan, and M. T. Rekveldt, Phys. Rev. A 69, 043615 (2004). 

[76] J. Erhart, S. Sponar, G. Sulyok, G. Badurek, M. Ozawa, and Y. Hasegawa, Nat. Phys. 8, 185 (2012). 



